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mixed-diffusion bistable equations 

Denis Bonheure, Juraj Foldes & Alberto Saldana*^ 


Abstract 

We consider a fourth-order extension of the Allen-Cahn model with 
mixed-diffusion and Navier boundary conditions. Using variational and 
bifurcation methods, we prove results on existence, uniqueness, positivity, 
stability, a priori estimates, and symmetry of solutions. As an application, 
we construct a nontrivial bounded saddle solution in the plane. 

Keywords: higher-order equations, bilaplacian, extended Fisher-Kolmogorov equation 

2010 Mathematics Subject Classification: 35J91 (primary), 35G30, 35B06, 35B32, 35B45 

(secondary) 


1 Introduction 


We study the following fourth-order equation with Navier boundary conditions 

( 1 . 1 ) 


— (3Au = u — 
u = Au = 0 


in n, 
on dQ, 


where V, C V > 1 is a bounded domain and /3 > 0. Such boundary condi¬ 
tions are relevant in many physical contexts [23) and they permit to rewrite (1.11 


as a second order elliptic system with Dirichlet boundary conditions. In our best 
knowledge (1.11 was analyzed only for V = 1, see [35] and references therein. 


In this paper we present results on existence, uniqueness, positivity, stability, a 
priori estimates, regularity, and symmetries of solutions in higher-dimensional 
domains when /? > -s/S- The case /3 < -s/S requires different approaches and 
techniques and we only prove partial results in this case. 


The problem (1.1) is a stationary version of 


dtu + — Au = u — u^ in 11 x [0, oo), 7 > 0, 


( 1 . 2 ) 
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which was first proposed in 1988 by Dee and van Saarloos as a higher-order 
model for physical, chemical, and biological systems. The right-hand side of 
(1.2) is of bistable type, meaning it has two constant stable states u = ±1 
separated by a third unstable state u = 0, see m The distinctive feature of 
this model is that the structure of equilibria is richer compared to its second 
order counterpart 


dtu — Au = u — in X [0, oo). 


(1.3) 


giving rise to more complicated patterns and dynamics. The equation (1.3) is 
related to the Fisher-KPP equation [Fisher-Kolmogorov-Petrovskii-Piscunov or 
sometimes simply called Fisher-Kolmogorov equationproposed by Fisher [20] 
to model the spreading of an advantageous gene and mathematically analyzed 
by Kolmogorov, Petrovskii, and Piscunov m- 


The equilibria of (1.3) satisfy the well-known Allen-Cahn or real Ginzburg- 
Landau equation 

— Au = u — in D (1-4) 

with the associated energy functional 




1 

' dx - 


— 1)^ dx, 


G H\n), 


(1.5) 


where Ff”(f2) = VF”’^(D) denotes the usual Sobolev space. This functional is 
used to describe the pattern and the separation of the (stable) phases ±1 of a 
material within the van der Waals-Cahn-Hilliard gradient theory of phase sepa¬ 
ration |12j . For instance, it has important physical applications in the study of 
interfaces in both gases and solids, e.g. for binary metallic alloys [T] or bi-phase 
separation in fluids [43]. In these models the function u describes the point- 
wise state of the material or the fluid. The constant equilibria corresponding 
to the global minimum points ±1 of the potential j(|mP — 1)^ are called the 
pure phases, whereas other configurations u represent mixed states, and orbits 
connecting ±1 describe phase transitions. 

To understand the formation of more complex patterns in layering phe¬ 
nomena—observed for instance in concentrated soap solutions or metallic al¬ 
loys—some nonlinear models for materials include second order derivatives in 
the energy functional. The basic model can be seen as an extension of (1.5), 
namely 

f [(V^it)^ -p 5 (m)|Vm|2 -p Win)] dx, u G H'^{n), 


where denotes the Hessian matrix of u. It appears as a simplihcation 
of a nonlocal model [28] analyzed in dimension one in [3 m ESI E7] and in 
higher dimensions in [I3l|22l|27]. In [27], the Hessian is replaced by Au as 
a simplification of the model and it was also proposed as model for phase-field 
theory of edges in anisotropic crystals in [33]. Finally, we also mention the study 

^in the original model, the nonlinearity it® is replaced by v? 
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of amphiphilic films [34] and the description of the phase separations in ternary 
mixtures containing oil, water, and amphiphile, see |25j . where the scalar order 
parameter u is related to the local difference of concentrations of water and oil. 


These models motivate the study of the stationary solutions of (1.2). After 
scaling, equilibria of (1.2) in the whole space solve 


A^u — PAu = u — in . 


( 1 . 6 ) 


We refer to (1.6) as the Extended-Fisher-Kolmogorov equation (EFK) for j5 > 0 
and as the Swift-Hohenberg equation for /3 < 0. This fourth-order model has 
been mostly investigated for N = 1, i.e.. 


— j3u' -I- — u = 0 


(1.7) 


When /3 S [-/S, oo), there is a full classification of bounded solutions of (1.7), 
which mirrors that of the second order equation. Specifically, each bounded 
solution is either constant, a unique kink (up to translations and reflection), or 
a periodic solution indexed by the first integral, whereas there are no pulses. 

For /3 S [0, ■\/8), infinitely many kinks, pulses, and chaotic solutions appear. 
Some solutions can be characterized by homotopy classes, but a full classifica¬ 
tion is not available. The threshold is related to a change in stability of con¬ 
stant states u = ±1. The proof of these results rely on purely one-dimensional 
techniques, for instance, stretching arguments, phase space analysis, shooting 
methods, first integrals, etc. For more details on the one-dimensional FFK we 
refer to [SlISH] and the references therein. 

For TV > 2 let us mention [5], where the authors prove an analog of the 
Gibbon’s conjecture and some Liouville-type results. Other related results can 
be found in [13 US mi- We also mention that the differential operator from 


(1.6) appears in geometry in the study of the Paneitz-Brans on operator^ see for 
instance USE], and as a special case of some elliptic systems, see e.g. [if)] . 
In the context of Schrodinger equations in nonlinear optics, the fourth order 
operator in (1.6) is used to model a mixed dispersion, see e.g. [MIT]. 


To introduce our main results denote 


H ■= H^{n)nH^{n) 


( 1 . 8 ) 


the Sobolev space associated with Navier boundary conditions (see [13 for ^ 
survey on Navier and other boundary conditions) and let Jp : H ^ M. be the 
energy functional given by 


^ , f f\Au\^ \Vu 




for u € H. 


(1.9) 


Any critical point u of Jg is a weak solution of (1.1), that is, u satisfies 
/ AuAv -\- jSVuS/v + (u^ — u)v dx = 0 for all v G H. 

Jq 
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Moreover, we say that u is stable if 

Jp{u)[v,v]= j |Ar;p + + (3r 

Jn 


— l)v^ dx >0 for all V G H 


and strictly stable if the inequality is strict for any v ^ 0. 


Extracting qualitative information even for global minimizers of (1.9) is far 


from trivial, since many important tools used in second order problems are no 
longer available. For example, one cannot use arguments involving the positive 
part M"*" := max{u, 0 }, absolute value, or rearrangements of functions since 
they do not belong to in general. Furthermore, the validity of maximum 

principles (or, more generally, positivity preserving properties) is a delicate issue 
in fourth-order problems and does not hold in general. 

For the rest of the paper, Ai(r2) = Ai > 0 denotes the first Dirichlet eigen¬ 
value of —A in n and a hyperrectangle is a product of N bounded nonempty 
open intervals. 

The following is our main existence and uniqueness result, for the second- 
order counterpart, we refer to [5]. 

Theorem 1.1. Let jd > 0 and 12 C with N > 1 be a smooth bounded domain 
or a hyperrectangle. If + PXi > 1, then u = Q is the unique weak solution of 

( 0 . If 

Xl + pXi<l, (1.10) 

then 

1) there is e > 0 such that (O admits a positive classical solution for each 

l-A? 


P G {l3 — e, P), where P = ^ 


2) for each P > there is a positive classical solution u of ( 1 . 1 ) 

such that ||M||L°o(n) < < 2 '^ tn 12. 


3) for every P > ^/S there exists a unique positive solution u of (1.1). More¬ 
over, this solution is strictly stable and satisfies 


\m\L^ 


qn) < 1- 


The smoothness assumptions on 12 are needed for higher-order elliptic regu¬ 
larity results. We single out hyperrectangles to use them in the construction of 


saddle solutions and patterns. Indeed, by reflexion, positive solutions of (1.1) 


in regular polygons that tile the plane give rise to periodic planar patterns. 

The quantities involved in Theorem |l.l| 2) are of a technical nature. Observe 
that (1.10 1 holds for all big enough domains. As mentioned above, the threshold 
'/S is related to a change in the stability of constant states m = ±1. For P > y/S 
the states are saddle-node type whereas for P < -s/S they are saddle-focus type. 
Hence in the latter case due to oscillations one can not expect u being bounded 
by 1. Such oscillations around one can be proved for radial global minimizers 
arguing as in 0 proof of Theorem 6 ]. Intuitively, for P > 1/8 the Laplacian 
is the leading term and the equation inherits the dynamics of the second order 
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Allen-Cahn equation, while for j3 £ (0, v^) the bilaplacian increases its influence 
resulting in a much richer and complex set of solutions. We present numerical 
approximation^ of positive solutions using minimization techniques in Figure 
[T] below. 



Figure 1: Numerical approximation of the global minimizer of ( |1.9[ l for Q. = [0,50]^ with 
/3 = 0.1 (left) and /3 = 4 (right). The dotted lines represent the level set {u = 1}. 


Note that Theorem o 1) holds for any /3 > 0, but only for appropriate 
values of Ai. 

Theorem o follows directly from Theorem |5.2| and Theorem |11.1[ The 
proof is based on variational and bifurcation techniques. For the variational part 
(Theorem |5.2[ ), we minimize an auxiliary problem for which we can guarantee 
the sign and L°° bounds of global minimizers. Next, we prove that global 
minimizers of the auxiliary problem are solutions to (1.91. The uniqueness 
is proved using stability, maximum principles, and bifurcation from a simple 
eigenvalue (Theorem |ll.l[ ) . After the paper was accepted Tobias Weth suggested 
us an alternative proof for uniqueness, see Remark 1 6.5 [ 

We depict a numerical approximatiorj^of the bifurcation branch in Figure m 
This branch can be continued even for /3 < 0, see Section m for an example of 
such a branch and we refer again to [31] for a survey on (1.7) for (5 < 0. See also 
Remark 11.2 for a discussion on the explicit values of the bifurcation points. 

We now use the solution given by Theorem |l.l| to construct a saddle solution 
for (O. We call u a saddle solution if u ^ 0 and u{x, y)xy > 0 for all (cc, y) £ 
See Figure]^ below. 


Theorem 1.2. For /3 > y' ^ the problem A^u — j3Au = u — m has 
a saddle solution. 


We refer to miiii] for more information on saddle solutions for second order 
bistable equations. 

^Computed with FreeFem++ m and Mathematica 10.0, Wolfram Research Inc., 2014. 
^Computed with AUTO-07P [18) . 
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Figure 2: Numerical approximation of the bifurcation branch and some radial solutions. 
Here is a ball in of radius 240.483. 



Figure 3: Saddle solutions for 0 < \/8 (left) and 0 > \/8 (right). The dotted lines are the 
level set {u = 1}. 


In the following we explore properties of positive classical solutions with a 
special focus on stability and symmetry properties. 

Theorem 1.3. Let LI be a ball or an annulus and let u be a stable solution of 
<0) with P > vT^ — 2Ai such that ||u||Loo(n) ^ 1- Then u is a radial function. 

Note that Theorem |1.3| does not assume positivity of solutions. We believe 
that the restriction on /3 is of a technical nature, but it is needed in our approach, 
see also Remark FTTI 

More generally, for reflectionally symmetric domains we have the following. 
We say that a domain is convex and symmetric in the ei — direction if for every 
X = {xi,X 2 , • ■ •, xn) G LI we have {(txi,X 2 , • ■ •, xn) '■ t G [—1,1]} C LI. 

Proposition 1.4. Let P > VS and let Ll C be a hyperrectangle or a bounded 
smooth domain which is convex and symmetric in the ei —direction. Then, any 
positive solution of O) satisfies 

u{xi,X 2 , • ■ ■, xn) = u{—xi,X 2 , ■ ■ •, xn) for all X = {xi,... ,xn) G Ll, 
dxiu{x) < 0 for all X G Ll such that xi > 0. 
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The proof of Proposition 1.4 follows a moving-plane argument for systems, 
see Figure (right) for the described symmetry. Note that for f3 S (0, -x/S) the 
solution oscillates when close to 1 in big enough domains, in particular it is not 
monotone, although it may still be symmetric. In balls, Proposition |1.4| implies 
Theorem|1.3|with the stability assumption replaced by positivity of the solution. 


In the following we focus on properties of particular solutions, namely, global 
minimizers. The next theorem states positivity of global radial minimizers, that 
is, functions u G Hr := {v G H : v is radial in fl} such that Jp{u) < Jp{v) for 
all V G Hr- 


Theorem 1.5. Let fl be a ball or an annulus, /? > 0, (I.IO) hold, and let u G Hr 
be a global radial minimizer of (1.9) with ||u||Loo(a) < I. Then drU does not 
change sign if LI is a ball and drU changes sign exactly once if LI is an annulus. 


Note that global minimizers satisfy the required bound for /? > \/8, see 
Proposition |4.I| below. The proof of Theorem |I.5| relies on a new flipping 
technique that preserves differentiability while diminishing the energy and it is 
therefore well suited for variational fourth-order problems. Theorem |I.5| clearly 
implies that global radial minimizers do not change sign and we conjecture that 
this property holds in general, even for (3 < i/S- For /3 large we can relax the 
assumptions on the solution, as stated in the following. 


Corollary 1.6. Let Ll be a ball or an annulus, ft > -s/l^— 2Ai, (1.10) hold, and 
let u be a global minimizer of (1.9) in H. Then u is radial and does not change 
sign in Ll. Moreover, drU does not change sign if Ll is a ball while drU changes 
sign exactly once if Ll is an annulus. 


Proof. From (3 > y/l2 — 2Ai and (1.10) follows /3 > and from Proposition 
|4.1[ we deduce that ||M||L°°(n) < 1- Theorem 
radial. Then u is in particular the global minimizer in Hr 
follows from Theorem 11.51 

□ 


therefore implies that u is 
and the corollary 


As last theorem we present a uniqueness and continuity result for 

'jA'^u — Au = u — in n, 

u — Au = 0 on dLl , 


( 1 . 11 ) 


when 7 —> 0 . 

Theorem 1.7. Let 12 C be a smooth bounded domain with the first Dirichlet 
eigenvalue Ai < 1. Let 7 > 0 and u.y be a global minimizer in H of 



iVup u'^ 

2 T “ yj 


dx. 


( 1 . 12 ) 


There is (5(12) = (5 > 0 such that, for all 7 € [0,(5], u.y is the unique global 
minimizer in H and u~^ > 0 in Ll. Moreover, the function <5] —>■ C"^(12); 
'y i-G' Uj is continuous and uq is the global minimizer of (1.12) in Hq(LI) with 
7 = 0. 
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The proof of this singular perturbation result relies on regularity estimates 
which are independent of the parameter 7 , see Lemma |2.3| 

The paper is organized as follows. In Sections Hi and 1 ^ we prove crucial 
auxiliary lemmas for obtaining a priori estimates of solutions. Sectionj^contains 
the proof of the variational part of Theorem 0 Our study of the stability of 
solutions is contained in Section and In Section we present our flipping 
method and prove Theorem |1.5| The proof of Proposition I.4| can be found in 
Section m and the saddle solution is constructed in Section l!?]! The bifurcation 
result involved in the proof of Theorem |1.1| is contained in Section [D and the 
continuity result, Theorem |1.7[ is proved in Section Finally, in Section [T^ 
we include two numerical approximations of bifurcation branches. 

To close this Introduction, let us mention that the case P < -s/S requires a 
different approach due to possible oscillations around 1. In particular, in this 
case strict stability, uniqueness, and symmetry properties of positive solutions 
and global minimizers are not known. 

Acknowledgements: We thank Guido Sweers and Pavol Quitter for very help¬ 
ful discussions and suggestions related to the paper. We also want to thank 
Christophe Troestler for introducing us to FreeFem-|--l-. 


2 Auxiliary lemma 


We prove a very helpful lemma that allows us to obtain a priori bounds on 
classical solutions. These arguments were used simultaneously in [6l Lemma 
3.1.] to obtain similar a priori estimates in . 

For the rest of the paper we denote by (70(11) the space of continuous func¬ 
tions in Q vanishing on dVl. 

Lemma 2.1. Let Lt C be a bounded domain, /3>0, / :M— J-Ka continuous 
function satisfying /(O) = 0, and let u G n Co(^) be a solution of — 

PAu = f{u) in 11 such that Au G C'o(^). Set u := rnaxu, u := rmnw, and 

o n 

g : K —>■ K given by g{s) := ^/(s) + s. Then 

u<maxg and naming. (2.1) 

[n,iA] 

Moreover, 

1. Ifu < maxq and f < 0 in (l,oo), then u < maxq. 

“ [o,H] V ) y 


2. Ifu> min a and f > 0 in (— 00 ,— 1 ), then u > min g. 

[«, 0 ] [- 1 , 0 ] 


Proof. Let w G (7^(11) fl (7o(ll) be given by w{x) := 


-Au{x) + fu(x). 


We 


prove only the second inequality in (2.1) as the first one follows similarly. Fix 
Xo,^o G H such that u(xo) = u and w(^o) = rmnui. If xq G 511, then 


; = u(xo) = 0 = g(u(xo)) > ming 




and the second inequality in (2.1) follows. If € dQ, then —Au + > 0 in 

n and the maximum principle implies that m = 0 and the second inequality in 
( 2 . 1 ) follows. Ifxoi^o G then w(^o) < w{xo) = — Au{xq)-\-^ u{x{)) < ^u{xo). 
Since —Aw{^q) < 0, 


/(w(Co)) 


-m(Co) = A^'u(^o) - ^Am(^o) 


-m(Co) 


= -Aw(^o) + < ^u{xo), 


that implies u > min g. 

[u-,u] 

We now prove claim 2. only as claim 1. is analogous. Assume m < — 1, 
otherwise the statement is trivial. Since / > 0 in (—oo,—1), then g(s) > s 
in (— 00 ,- 1 ), and therefore min q > min s = u. Thus, if u > mino then 

[u.-l] ^ ^ [ti.O] 

u > min q as claimed. 

^ [- 1 . 0 ] 

□ 


2.1 Regularity 


In this section we prove two regularity results. The first one is a rather standard 
application of known arguments in the fourth order setting and we include 
details for reader’s convenience. The second lemma is more subtle and the 
crucial point is the dependencies of the constants involved. 

Recall the definition of the space H given in (1.8). 


Lemma 2.2. Let 12 C N > 1 be a smooth bounded domain or hyperrect¬ 
angle and fix 13, "f > 0, and f G Let u G H be a weak solution of 

■jA'^u — /3Au = f, that is, 


-fAuA(t) +(3\7u\7(l)-f(l)dx = 0 'icjiGH. (2.2) 

Jn 

Then for each p > 1 one has u G H with u = Au = 0 on dLl and 

l|ii||iV'‘.p(n) < C'll/l|L°°(n). 

where C = In addition, if f G for some a G (0,1), then 

u G C'4’“(n) with Au G Coin) and 

lhllc«.“(a) < (5||/llc“(n)) (2-3) 


where C = C{(3,^,p,Tl). 

Proof. Throughout the proof C > 0 denotes possibly different constants de¬ 
pending only on (3, 7 , p, and 12. 
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Assume first that SI C N > 1, is a smooth bounded domain. By 

the Riesz representation theorem there are weak solutions u,v G i?g(Sd) of the 
equations 

—Av + —V = / and — 'jAu = v in SI. (2-4) 

7 

Fix p > 1. Then, by [24l Lemma 9.17] and [33l Ch 9 Sec 2 Thm 3], 
||M||vF'‘>p(n) < C'||il||rF 2 ,p(f 2 ) < C'||/||Lp(n)j which gives the first estimate in the 
statement. 

By embedding theorems, we have u G C^(J1) and u = 0 = v = —7 Am on 
dfi. Let M G iL be as in the statement. By integration by parts, u satisfies \2.2\ 
and M = u in S7, since weak solutions of \2.2\ are unique. Finally, if / G (^“(Sl), 
by Schauder estimates |32l Thm. 6.3.2], v G u G for some 


a G (0,1), and ( |2.3| ) holds. 

Now without loss of generality assume that is a hyperrectangle of the 


form n = nLio, k] for some k > 0, i = 1,..., A. Let u,v G Aq (fl) be weak 
solutions of ^2.4] ). Using odd reflections and the Dirichlet boundary conditions 

tN 


we extend u and v to Cl = Y\i^i[—li,2.li] and obtain weak solutions of (2.41 
(with the odd extension of / as well) defined in VL. Then interior regularity [331 
Theorem 1 Sec 4 Ch 9] implies that for any fli C 142 C 

l|u||iv^.p(ni) < C'(||/||L“(n 2 ) + 11^11^2(02)) • (2-5) 

Note that we replaced by L? on the right hand side, which can be done 


using Sobolev embeddings and iteration in (2.5). Also, ||u||i 2 (Q) ^ ^II/IIl 2 (q) 
by testing (2.4) by v and u respectively and using standard estimates. 


As before, m = it in 14 by integration by parts and uniqueness of weak 


solutions. Finally, (2.3) follows analogously from interior Schauder estimates 
[321 Thm. 7.11]. 

□ 

Lemma 2.3. Let LI be a smooth domain, / 3,7 > 0, and let u G H D L°°{L1) 
be a weak solution of — fdAu = u — in H. For every p > 1 there is 

7 o = 7 o(/ 3 )P)^ 2 ) > 0 and C = C{P,p,Ll) such that i /7 G (0,7o) then 


|w||rve-p(n) < <^(1 + lhlli~(o))- 


( 2 . 6 ) 


Proof. In this proof, C denotes different positive constants which depend on 
14,/3, and p, but are independent of 7 . 

that u G (7^’“ (14) 


Since u G L^{Ll) we have by bootstrap and Lemma 


2.2 


with Au G C‘^(L2). Then (u,v) with v = —jAu and f = u — solves (2.4) in 
the classical sense. By [331 Ch. 8 . Sec. 5. Thm 6 ], there is 70 = 7 o(/ 3 ,p, 14) > 0 
such that, for every 7 G (0,7o), 

Ikllw2.p(n) < <^(1 + ||u|||c=o(n)). 

On the other hand, —/3Au = u — — yA^rt = u — + Av and by [Ml Lemma 

9.17] for every p G (0, 00 ) we have 

I|w||iv 2 .p(n) < C\\u — + Am||lp(o) < C(1 + ||M||ioo(Q)). (2.7) 
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In particular, ||AM|| 2 ,P(n) < C(1 + Set w = Au. Since u G 

we have that re is a weak solution of 


'jA'^w — f3Aw = fi := A(u — u^) = w — 6it|Vitp — 3u^w in fl 


with w = 0 on dfl. Moreover, Aw = A^u = ^{u — + fiAu) = 0 on 917, since 

u G C^’“(r2) and u = Au = 0. Observe that ||/i||Lp(n) < <^(1 + IKIIloo(q))- 
Then we can repeat the argument for w and /i instead of u and / respectively 
to obtain that w G C^'°'(Vl) and 

Il^i'lliv2.p(n) < C'Cl + lkllL“(n)) ■ (2-8) 


Then ^ := Aw G (^^’“(O), ||C||LP(n) < 0(1 + and ^ is a weak solution 

of — j3A^ = /2 := A/i and ^ = 0 on 9f2. Additionally, A^ = A^w = 
^(/i + 13Aw) = 0 on 90, where we used u = Au = w = Aw on 90. Here 
we are fundamentally using that fi vanishes on 90. Note that ||/ 2 ||LP(n) < 
(7(1 + ||n||® oo(f 2 ))- Thus we can iterate the above procedure one more time with 
u and / replaced by ^ and /2 respectively, to obtain 


IICIIiv2.p(n) < (7(1 + ||u||i<»(n)) • 


(2.9) 


Note that we cannot iterate anymore since /2 is not vanishing on 90, and 
therefore we cannot obtain boundary conditions for higher derivatives. Finally, 
(Tel follows by (ITtI), d^, d^, and El Ch.9 Sec. 2 Thm 3]. 

□ 


3 A priori bounds 

For any /3 > 0 we fix the following constants for the rest of the paper 

'4 + /32 ' ' 


Cp:= 


4 + /32 




Kn := 


^ 27-2 


1.58. (3.1) 


Remark 3.1. IFe now state some facts that will be useful later in the proofs, 
i) Cp > 1 is the unique positive root of h{s) = -^{s — s^) + s (cf Lemma 


2.1). In particular, h > 0 in (0, (7/3). 


ii) Mp > 1 is the maximum value of h in (0, oo) 

Hi) If f3 > 's/S then h is increasing on (0,1). 
iv) Mp <Cp if /3> Kq. 

Lemma 3.2. Let O C be a bounded domain, j3 > 0, and let u he a classical 
solution of (1.1). 

i) If u is nonnegative in O, then ||u||Loo(n) < Mp. 
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ii) If 13 > and u is nonnegative in O, then < 1 . 

Hi) If (3 > Vs and ||u||ioo(n) < \j^ + 3, then ||M||Loo(f 2 ) < 1. 

Proof. We use the notation of Lemma 2.1 with f{s)=s — s^. We prove claim 
Hi) first. Assume without loss of generality that ||M||L=o(f 2 ) =u< + 1. By 

direct computation g{s) < —s for s S (—^ + 1, 0). Therefore 


max g(s) < max —s = u. 

sG[-u,0] sG[-«,0] 


(3.2) 


2.1 


u < max g < max g. Then, in virtue of (3.2), we get 
from Lemma [ 2 Tl and Remark |3.1| iii) that u < max 5 = 1 


i |2.1| and Remark |3.1|iii) 


[ 0 , 1 ] 


Moreover, by Lemma 

u < max q. It follows 

[0,H] 

as claimed. 

Now, let /3 > 0 and assume u > 0. By Lemma |2.1| and Remark |3.1[ u < 
max g = Mg and the first claim follows. If /3 > v^, then g is nondecreasing in 

[O.oo] 

[0,1] and, by Lemma |2.1l u < max g = 5 ( 1 ) = 1. 

□ 


4 Bounds for the global minimizer 


Recall the definition of ATo, Mg, and Cg given in (3.1) and of H given in (1.8). 


Proposition 4.1. Let v be a global minimizer of (1.1) in H. Then 

\\v\\L’=°{n) < Cg ifl3>Ko and ||u||Loo(n) < 1 if P > Vs. 

When (3 G (0, -s/S) we expect that global minimizers are not bounded by 1 
in big enough domains, see Figure and Figure 


Proof. Assume /3 > Kq, let Jg be as in ( |1.9[ ), and let v denote a global minimizer 
in H, i.e. Jg{v) < Jp{w) for all w G H. We show first that v G Let 

1 sign(s)(C '/3 — Cg) = —^ sign(s)C '/3 otherwise , 


f-.R- 
and let 


J :H 


J{u) := f 
Jn 


|Ar 


V^^-Fiu)dx, (4.1) 


where F{s) := fit) dt. Let u be the global minimizer of J in H. By Lemma 
12.21 u is a classical solution of 


— pAu = f{u) in LI, 

u = Au = 0 on dLl. 


(4.2) 
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Using Lemma|2.1| (and its notation), since g > 0 in (0, oo), we have u > min g = 

I_I [m.U 

ming, and consequently u > min g > —Mg. If/3> we additionally have 
[«. 0 ] ^ [- 1 , 0 ] 
u > min g = — 1. 

“ [- 1 , 0 ] 

Replacing u by —u and noting that 13 > Kq we conclude ||M||Loo(n) < Mg < 
Cg and ||u||ic=o(Q) < 1 if /3 > -s/S- In particular f{u{x)) = u{x) — u^{x) for 
X G n, and therefore 


J{u) = Jg{u) > Jg{v) > J{v), 

where the last inequality is strict if 111111 ^= 0 ( 0 ) > Cg, a contradiction to the 
minimality of u. Thus v G and, by Lemma 2.2 u is a classical solution 

of (4.2 1 . If /3 > -x/S and ||u||ioo(Q) > 1 , then, by Lemma 3.2 part in), ||u||ioo(Q) > 


\J^ 3- 1 > Cg and we obtain a contradiction as above. 


□ 


5 Existence of positive solutions 


Recall the definition of Kq, Mg, and Cg given in (3.1). 
Lemma 5.1. Let U C he a bounded domain, /3 > 0, 


/: 


f(s) 

Cg-Cl = -^Cg 


ifs<0, 
if s G [0, Cg], 
if s> Cg, 


(5.1) 


and let u be a classical solution of (4.2 1 for this choice of f. Then 0 < u < Mg. 
Moreover, if (3 > y/S, then u < 1. In particular, if j3 > Kq, then u is a classical 
solution of ( 1 . 1 ). 


na |2.1| (i 
, that IS, 


(and using the same notation) u > min g > 0 by the 


u > 0 in U. On the other hand, again by Lemma 2.1 


Proof. By Lemma 

definition of C,g 
we have that u < max g = max q = Mg. 

]0.u] [0.1] 

If /? > -s/S, then g is nondecreasing in [0,1], and therefore u < maxg implies 


by Lemma 


2.1 


[o,«] 

that u < max g = g(l) = 1. Finally, if /3 > Kq, then Mg < Cg, 


and thus f (u(x)) = u{x) — u{x)^ for all x G that is, u solves 


□ 


Theorem 5.2. Let fi > 0 and U C with V > 1 &e o smooth bounded domain 
or a hyperrectangle. If + /3Ai > 1, then u = t) is the unique weak solution of 
(1.1). If Af + l3Xi < 1, then for j3 > Kq there is a positive classical solution u 
of (1.1) such that ||it||L“(n) ^ Mg and Au < in U. Additionally, if [3 > \f8 
then ||u||L=o(n) < 1 . 
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Proof. Let + /3Ai > 1 and assume by contradiction that there is a nontrivial 
weak solution w G iL of (1.1). Testing equation with u yields 


0= / IAm|^ + ^|Vu|^ + dx > (A^ +/3Ai — 1) j u^dx>0, 

Jq. j Vt 

by the Poincare inequality, a contradiction. 

Now, assume that Xf + /3Ai < 1, let / be as in (5.1), and J, F as in (|4.1| 

for all s G 


|n 


for 


with this choice of /. Note that F(s) < | for all s G K. Thus J(u) > — 4 
all u G H. Standard arguments show that J attains a global minimizer u in H 
and that u is a weak solution of A^u — /3Au = f{u) in with Navier boundary 
conditions. Observe that m ^ 0, since Xf + /3Ai < 1 implies J(i5<pi) < 0 for 


sufficiently small S 
Laplacian in fi. 


> 0, where is the first Dirichlet eigenfunction of the 


Arguing as in Proposition 4.1 we have t hat u G and therefore u G 

Then u is a classical solution 


Pi (7(11) and Am G (7o(11), by Lemma 


2.2 


of (1.1) satisfying 0 < m < Mp if ,5 > Kq, and 0<M<lif/3> -s/S, by Lemma 

EH 


The strict positivity of u (recall u ^ 0) and — Au+ is a consequence of the 
maximum principle for second order equations and the following decomposition 
into a second order system 


A,/3 3 

—Au+—u = w, — Am;+— ui = (lH ——)u — u 
2 2 4 


in n, u = w = 0 on 912, 


where (1 + ^)u — > 0 since 0 < w < Mp < Cp for /3 > Kq. 


□ 


6 Stability of positive solutions 


In this section we prove the following result. 

Theorem 6.1. Let dfl be of class ( 7 ^A ^ Then any positive solution 


of ( 1 . 1 ) is strictly stable. 


The proof of Theorem 6.1 is an easy consequence of the following. 

Proposition 6 . 2 . Assume that ft > -s/S, o.nd let u be a positive solution of 
( |1.1[ ). Then, 

L I Aup + /3|Vup + (m^ — l)z;^ dx 
hi = inf - r —- = 0 . 

veH\{o} j^v^dx 

Indeed, assume for a moment that Proposition |6.2| holds. Then we have. 
Theorem \6.1[ Let TL := {v G H : ||M|p 2 (Q) = I}. By Proposition 


6.2 


inf [ |Avp +/3 |VmP + (3u^ — 1 )m^ dx = + 2 inf [ dx > 0 


and Theorem 16.11 follows. 


□ 
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The proof of Proposition |6.2| uses the following result adjusted to our situa¬ 
tion. 


-A 0 


0 -A 


and M be a 2 x 2 continuous matrix 


Theorem 6.3. Let L := 

such that 

1) —M is essentially positive, that is, —M 12 > 0 and —M 21 >0 in LI. 

2) M is fully coupled, that is, Mi ^2 ^ 0 and M 24 ^0 in LI. 

3) there is a positive strict supersolution (p of L + M, i.e., a function (f > 0 
such that {L -|- M)(j) > 0 in LI. 

Then there are v,w G p| C'o(fl) unique (up to normalization) posi¬ 

tive functions such that 


(L + M) ( 1) = Ao , 

\w / \w 


( 6 . 1 ) 


where Aq > 0 is the smallest eigenvalue (smallest real part) of ( 6 . 1 ). 


Moreover, there are positive functions v,w G fj (70(11) unique up 

to normalization such that 


iL + M)[l)=X,Br 


( 6 . 2 ) 


where As > 0 is the smallest eigenvalue (smallest real part) of (6.2) and B ^ 0 
is a matrix with Bij G C{Ll) and Bij > 0. 


Theorem 6.3 is a particular case of [m Theorem 1.1] and [3S1 Theorem 5.1]. 
In [Sni Theorem 5.1] the result is formulated for matrices B G Co(L2), but the 
same proof applies for B G C(L2). 


Proposition \ 6. 2\ Let 


L := 

and 


-A 0 
0 -A 


M := 






B := 


0 0 
1 0 


, (6.3) 


~ 52 / /3 _i 


Let (/?! > 0 be the first Dirichlet eigenfunction of the Laplacian in LI, i.e., 

—Aipi = Xiipi in LI, (fii = 0 on dLl. (6-4) 

Note that —M is fully coupled, M is essentially positive by Lemma 3.2 ii), and, 
for /? > -s/S, the function (p = ^ positive strict supersolution of L -I- M, 
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because Ai + f — 1 > 0 and Ai + — 1 + f > 0. Then, by Theorem 


6.3 


there 


are unique (up to normalization) positive functions v,w G p| C'o(n)) 

such that 


(L + M) 


= {Xb-^)B 


where \b > 0 is the smallest eigenvalue. 

Since dVl is standard regularity arguments imply that v G C"^(n) n 
C'o(^), An G Co(ri) and then 

fl2 

A^n —/3An + — l)n = (As- —)v = y.iv in n, (6.5) 

^2 ■ ■ 

where ni := Xb — ^ is the smallest eigenvalue of (6.5). By multiplying this 
equation by u > 0 and integrating by parts we get that 


0 = ^1 


; dx 


In 


since u is a solution of (|1.1|). Therefore /ii = 0, because u and v are positive. 

□ 


This ends the proof by the variational characterization of fj,i. 


Remark 6.4. Note that without the assumption 0 > the result still holds for 
any solution u such that 0 < u < 1 inil, Ai + ^ — 1 >0, and Ai — 1 + ^ > 0. 

Remark 6.5. Uniqueness of positive solutions of (1.1) when 0 > -x/S can he 
proved using Proposition^ 


_ \as follows. Let u and v denote two positive solutions 

of ( |1.1[ ). Then w := u — v solves A^in — 0Aw + (u^ — l)w = —(uv + v'^)w in fl 
and w = Aw = 0 on dfl. By testing with w we have 

/ \Aw\'^ + 0\\/w\'^-\-(u^ — l)w'^ dx = / (—uv — v^)w‘^dx<0. 

Jn Jn 

But then Proposition\6.S\ yields that w = 0. 


7 Radial symmetry of stable solutions 


Theorem ES Let L and B be as in (|6.3|) and let 


Q:= 




Note that 


Q '■— Q + — 2^ B — 
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is a fully coupled matrix and —Q is essentially positive (see Theorem 6.3). 
Moreover, since (3 > \/vl — 2Ai there is i5 > 0 such that l/(Ai + ^) < <5 < 
(Ai + f )/3. For such 5, if is as in (6.4), we have 


(L + Q) 


<Pi 


((5Ai + - l)pi 

(Ai + 3(u^ — 1)(5 + f)(/?! 


> 0 , 


and therefore L + Q has a positive strict supersolution. Therefore, by Theorem 


6.3 


there is As > 0 and positive functions v,w G (il) H (70(^2) such that 


(L + Q) 


w 


= p B 


(3^ 


with pL = \b -^+2. 


By standard regularity arguments, v G H C'o(n) and v solves 

— /3Au + (3 m^ — 1 )?; = pv in Q, 

that is, V is the first eigenfunction and /r is a simple (first) eigenvalue. Stability 
of u implies /i > 0. We now argue by contradiction, assume that u is not 
radially symmetric. Then there is a nontrivial angular derivative ue = dgu ^ 0. 
In particular, ug must change sign in and A^ug — PAug + (3u^ — l)ug = 0 
in n, since dg and —A commute and u is a solution of (1.1). This implies that 


Ug is a sign-changing eigenfunction associated to the zero eigenvalue, but this 
contradicts the fact that ^ > 0 is the simple first eigenvalue. Therefore u is a 
radial function. 

□ 


Remark 7.1. Note that Theorem 1.3 holds true for any Ai if P > \/T2. On the 
other hand, if u is a nontrivial solution, then, by Theorem \5.^ we have that 
Xi + pXi < 1. Combined with P > '/l2 — 2Ai, we obtain that the infimum of 
P’s satisfying these inequalities is -s/S with corresponding Ai = -s/S — 


8 Positivity of global radial minimizers 


Before we prove Theorem |1. 5 [ we introduce some notation. Let 


r [0,R) iin = BR, 

\ {Ro,R) iin = BR\BR„ 


for some R > Rq >0, where Br denotes the open ball of radius r centered at 
the origin. To simplify the presentation we abuse a little bit the notation and 
also denote u{r) = u{\x\). Recall that Hr = {u G H : u is radially symmetric}. 

Theorem m Suppose first that fl = Br{0) and for a contradiction, assume 
that u changes sign. Then either u or —u has a positive local maximum in 
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(0,i?). Without loss of generality, assume there is r] G (0, i?) such that 1 > 

M := uiri) = maxu > 0. Let w : / —?> M be given by 

^" [0,K] 


V := 


\j^{M-u)+M in[ 0 ,? 7 ], 

u in (rj, i?] . 


Note that v is just a rescaled reflection of u with respect to u = M in [0, ry]. 
Since ^'(ry) = u'(ry) = u'(0) = u'(0) = 0 we have that v G C^(/). Also 

\v'\ < \u'\ in (0,i?) and \v''\ < \u"\ in (0, i?) \ {ry}, 
with strict inequalities in (0,ry), whenever u' ^ 0 and u" ^ 0. 


Furthermore in (0, ry) one has 0 < u < {M + 1) + M = 1 and 
'1 - M 


1 + M 

l-M 


+ 1 ) (M — u) > 0, 
l-M 




1 + M 


+ 1 M = 


2M 
1 + M 


(u + 1) > 0. 


Then — v? = {v — u)(v + u) > 0, thus |1 — = 1 — < 1 — = |1 — 

in [0, i?], and then 


J- 1 )^ dx < j{u^ - 1 ) 


— 1 )^ dx. 


( 8 . 2 ) 


By (8.11 and (8.2) one has Jp{v) < Jp{u), a contradiction to the minimality of 
u and thus u does not change sign in Q. 

The proof for the annulus H (0) for some R > Rq > 0 is 

similar. If u changes sign, then there are t], fj, G I such that 


M := u(r]) = maxu > 0 and m := —u(n) = — minu > 0. 

10,R] [0,i?] 


Without loss of generality, assume that ry < yr. Let u : J —>■ M be given by 

u in [i?o,?7]) 

v.= i ^^{M-u)+M in(ry,/i). 


—u 


Then in (?y, /i) 

m — M 


V < 


V = 


V + u = 


m + M 
M — m 
m + M 
M — m 
m + M 
M — m 
m + M 


in [fj., i?]. 

(M + m) + M = m<l ifTO> M, 


{u — M) + M<M<1 ifm< M, 


- Ij u + 
+ 1 I M + 


m 


- M 


m + M 
m — M 
m + M 


+ 1 M = 


2 m 


+ 1 M = 


m + M 
2M 

m + M 


(—u + M) > 0, 
(u + m) > 0 . 
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Therefore (8.1) (with strict inequality on and (8.2) also hold in this case. 

Then Jp{v) < Jp{u), which contradicts the minimality of u and thus u does not 
change sign in Q. 

It only remains to prove the assertions about drU. We prove the case when 
n is an annulus, since the case when is a ball can be treated similarly. 

Without loss of generality assume m > 0 in fl. Now, suppose by contradiction 
that u' has more than one change of sign. Then u' has at least three changes 
of sign: two local maxima and one local minimum. Let ryi € / be such that 
u{rji) = maxM =: Mi and fj £ I another local maximum. We only prove 

the case rji < fj, the other one is analogous. Let ^ G {r]i,fj) be such that 
u(/r) = min u =: m and 772 G such that 77 ( 772 ) = max it =: M 2 . Clearly 

m < M 2 < Ml. Define t; : / —>■ [0,1] by 


u 

^■■={ %^{n-Mi) + Mi 

M 2 


in [i?o,?7i] U [ 772 , 1 ?], 
in (?7 i,m), 
in [^ 7 , 772 ), 


By similar calculations as above it is easy to see that v G C'^(J), and that 
Jpiv) < Jp{u), which contradicts the minimality of u and thus u' only changes 
sign once in I. 


□ 


9 Symmetry of positive solutions 

Proposition \l .4\ We can write ini> as the following system. 


A 

— AtI + —71 = W 


/32 


\ h' ^ M 

—Aw + —w = u — u -\——u 
2 4 

u = w = 0 


in D, 

in D, 
on dfl. 


(9.1) 


By Lemma 3.2 and (3 > one has ||77||ic=o(f2) < 1 and it is easy to check that a 
standard moving plane method as in [12] can be applied to (9.1). The regularity 


of the boundary can be relaxed up to Lipschitz regularity if one uses maximum 
principles for small domains, see mm- 

□ 


10 Saddle solution 


Recall the definition of Kq, 


Mfj, and Cp given in 


(3.1). 


Theorem El Let /3 > Ko. 


By odd reflection, it suffices to find a positive 
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u G C"^(K+) solving 


— /3Au = u — 
Au = u = 0 




on 




( 10 . 1 ) 


where = {{xi,X 2 ) G X\ > 0 ,a ;2 > 0}. Indeed, denote again by u the 
extension to by odd reflections. Since u = Q on Hi := {x : X 2 = 0} one 

= 0 


= 0 on i7i. Also Am = 0 on Hi implies 


= 0 on Hi. Since u. 


Xu Uxixij 


U, 


Ux 


has Uxi — Uxixi — '^x 
on i7i, and consequently Ux 
Ux 2 X 2 Xi are odd with respect to Hi they are continuous on Hi. Moreover, all 
other partial derivatives (up to third order) are continuous on Hi as well, since 
they are even functions. The same procedure applies to H 2 := {x ■. xi = 0} and 

thus u G C3(K2). 


From the equation (10.1) we also obtain continuity of A^m and A^m = 0 on 


_i_. By a similar reasoning as above we obtain that the extension is of class 
C'"*^(M^) and it is a classical solution of < 0 ) with n = and consequently it 
is a saddle solution. 

We find a positive solution of ( 10. 1[ ) by a limiting procedure using Theorem 
5.2 with := (0, RY and letting R —>■ 00 . Note that for R big enough, the 

solution Ufi given by Theorem 5.2 satisfies that 0 < ur < Mp in and 


therefore, by Lemma [2^ there is some C > 0 independent of R such that 


Ikflllc4.o(o„) < C for all i? > 0 . ( 10 . 2 ) 

By the Arzela-Ascoli theorem there is a sequence Rn —>■ 00 such that urj.^ —)■ 


u in C'^(IR^) for some u satisfying (10.1). We now prove that m ^ 0. Indeed, fix 

r > 80(1 + ^)C'2 + 160, (10.3) 

where C > 0 is a constant independent of R specified below. We show that 

l|Mi?||L=o([0,r+2]2) - ^ for all i? > r + 3 . 

Assume by contradiction that there is i? > r + 3 such that 

1 

7^' 


ll'Ufl||L“([0.r+2]2) < 


(10.4) 


(10.5) 


We define the following sets 


Wi 

a;2 

W3 

1x4 


{x G 
{x G 
{x G 
{x G 


dist(a;, B^r) < 1 , xi < r + 1 , a ;2 < r + 1 }, 
dist(a:, B^Ir) > 1 , < r + 1 , a ;2 < r + 1 }, 

r + 1 < max{xi, 0 : 2 } < r + 2 }, 
a;i > r + 2 , or 0^2 > r + 2 }. 
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Note that Qr — lJi=i‘^*- Now, let (/>! e C'^(ni^) p] C'o(riK), (/)2 G (7^(0/?) such 
that 0 < (/)i < 1, ||^i>i||c2(nH) < itl for i = 1,2 and some K > 0 independent of 
R, and 


(pi = 1 in a; 2 , pi = 0 in W 4 , p 2 = 0 in W 2 , p 2 = ^ in a; 4 . 

Further, let p S P ( 70(11 i?) be given by tp := pi + p2UR. Then = 1 in 

102 , '0 — ^-R in W 4 , and there is some (7 > 0, depending only on C (from (10.2)) 
and K, such that ||0||c2(nij) < C. For i = 1,... ,4 let 


Ji{v) := f 

J U), 


+ jiv'^ - 1)^ dx for V e 


Note that J2 for ^ ^ here is as in (1.9l for 


2 = 1 


Q = Qft. Then J2 M^) < [(^ + + l](l^il + I'^sl) + J4{ur), and by (|10.5|), 

2=1 

I] Mur) > ^ + JMr). Therefore 


2=1 


Mur) - Mi’) + 2)) > 0, 


by (10.3), a contradiction to the minimality of ur. Therefore (10.4) holds and 
the maximum principle yields that m > 0 in is a solution of ( 10 . 1 ). 

□ 


11 Bifurcation from a simple eigenvalue 

Theorem 11.1. Let LI C K^, N > 1 be a smooth bounded domain or a hyper¬ 
rectangle. If the first Dirichlet eigenvalue Ai < 1, then there is e > 0 sueh that 
O) admits a positive solution up G C^’°‘(fl) for all /3 € (/3 — e,/3), where 

B = ( 11 . 1 ) 


Additionally, if P > VS and fl is smooth, then (1.1) admits a unique positive 
solution Up such that ||u, 8 ||L°°(n) < 1 for all /3 G [V 8 , P). 


Proof. Let X = {u G (7"‘’“(n) (7 (7^(S1) : u = Au = 0 on 911} and Y = (7°’“(H). 
Consider the operator 


(7 : K X X —>■ F; G{P, u) := A^u — PAu — u + u^. 

Then we have G{P, 0) = 0 for all p. Moreover u G X solves if and only if 
G{P, u) = 0. We consider the partial derivative 

duG : (0, 00 ) X X —>• £(X, Y), 9„(7(7,M)[n] = — pAv — v + hu^v. 
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For /3 > 0 we set Ap := 9„G'(/3,0), i.e., Apv = — /3Ar; — u. Let N{Ap) and 

R{Ap) denote the kernel and the range of Ap respectively. Note that v S N{Ap) 
if and only if A^v — /3Av = ?; in fl. Let tpi be the first eigenfunction of the 
Laplacian in fi, see ( |6.4[ ) with ||</ 5 i|| l 2 (q) = 1. 

By the definition of /3 > 0, one has ipi G N(Ap). Moreover, by the Krein- 
Rutman Theorem N(Ap) = {aipi : a G K}. Further, since Ap is self adjoint, 
by the Fredholm Theory R(Ap) = {v G Y : J^ipivdx = 0}. In particular, 
-^Apipi \p=p= —Aipi = Aitpi ^ R{Aj). Hence, by [13 Lemma 1.1] there are 
e > 0 and —functions fi : (—£,s) —)■ (0,oo) and u : {—£,£) -G X such that 
/3(0) = P and G{/3(t),u{t)) = 0 for all t G (—£,£), moreover G“^({0}) near 
(,0,0) consists precisely of the curves u = 0 and (/3(t), u(t)), t G (— e,e). Since 
duuG($,0)[(pi,(pi] = 0 and ^^(piduuuG{P,Q)[ipi,ipi,ipi]dx = 6 ipf dx > 0 we 
have a subcritical bifurcation, and therefore u = ±c(/3 — /?) 5 + o{t 2 ) for some 

constant c > 0 and t G (0, e). This proves the first claim. 

For the second claim, assume jd > -s/S and that H is a smooth domain. 
Let (0,r) be the maximal interval of existence in [-s/S, ,0] for the curve u with 
T G (0, ooj. By this we mean that /3(t) > -s/S for t G (0,T), that 7 (t) can be 
uniquely extended for each t G (0,T), and, if (3{T) > -s/S, then 7 cannot be 
uniquely extended at T. In particular, the curve ceases to exist if it intersects 
another curve e.g. (/3,0) or if it bifurcates. 

We show first that u(t) > 0 for all t G (0,T). By the G^— continuity of 
the curve C := {u{t) : t G (0,r)} we have that 0 < < 1 for all t G (0,T) 

sufficiently close to zero. By Lemma 3.2 iii) and the continuity of C it follows 
that ||n(t)||^=o(Q) < 1 for all t G (0,T). To show that u{t) > 0 for all t G 
(0,r) we argue by contradiction. Assume that i = sup{t G (0,T) : u(s) > 
0 in for all s G (0,t]} < T. Let u = u(t) G G^’“(H). Note that u satisfies the 
system 


—Au fiu = w in fl. 


—Aw = u — 


in 


w = u = 0 on dfl, 


for some /3 > 0, and u — > 0 since 0 < ft < 1. Since T > t we have that 

u ^ 0. Then, by the maximum principle and the Hopf Lemma, u > 0 in H and 
< 0 on dil, where ly denotes the exterior normal vector to dil. But this 
contradicts the G'*—continuity of C and the definition of u. Therefore u{t) > 0 
for all t G (0, T). 

Then, by Lemma 3.2 ||u(t)||L~(n) < 1 for all t G (0,T), and standard 
elliptic regularity theory yields that ||M(t)||c^'“(n) ^ C' for all t G {0,T) and 
for some G > 0. Moreover, since the positivity is preserved along the curve, 
u{T) ^ 0. Indeed, {/3{t),u{t)) cannot return to a neighborhood of (,0,0) by 
uniqueness close to (,0,0). Also, it cannot intersect (/?,0) as any other branch 
bifurcating from (,0,0) (/3 < ,0) consists locally of sign changing solutions because 
the corresponding eigenfunction directions are sign changing (perpendicular to 
the principal eigenfunction). By the first part of Theorem 5.2 we know that 
/3(t) < j3 for all t G (0,T). This implies that necessarily (-s/S, (3) C {,0(t) : t G 
(0,r)}. This proves the existence of solutions for all (3 G [-s/S, ,0). 

We now show that up is the unique positive solution of ([gfor/3G [V 8 ,^). 
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Indeed, let v denote a positive solution of (1.11 for some /3o & [a/S,/ 3). By 


Theorem 6.1 n is a strictly stable solution, and therefore DuG{v, /3o) is an 
invertible operator. Then, by the implicit function theorem there exists e > 0 
and a smooth curve 7 : (/3o — e,/3o + s) X such that = 0 and for 

any solution of G{/3,u) = 0 sufficiently close to {/SojV) one has m G 7. 

Arguing as before, we can extend 7 to a maximal interval {Pi,/32) with 7 
containing only positive solutions. Then the strict stability in Theorem |6.1| 
guarantees that 7 does not have bifurcation or turning poi nts. Since the only 
solution ioi f3 > ^ is zero, by the first part of Theorem 5.2 we have that 


132 < P and 7(/32) is a non-negative function. Arguing as before, one obtains 
that necessarily ^{(32) = 0 and then /32 = Here, as above, we have used that 
all other bifurcation points of the form (/3, 0) must correspond locally to sign 
changing solutions. The uniqueness of the branch close to the bifurcation point 
(,d,0) yields that necessarily v = as desired. 

If H is a hyperrectangle, one proves the positivity along the curve using 
Serrin’s boundary point Lemma |39[ Lemma 1] at corners and the rest of the 
proof remains unchanged. 

□ 


Remark 11.2. For balls of radius R > 0 we can explicitly write the relationship 

between R and the bifurcation point (3. Indeed, in this case, R := N/ 2 - 1,1 ^ 

y —P+\/ ^^+4 

where N is the dimension and jAr/2-1,1 the first positive zero of the Bessel 

function Jjv/2-ij see for instance \2tA Section f.l]. For example, for (5 = y/%, 

the bifurcation occurs at balls of radius Rn ■= for instance, R 2 

V V3-V2 

4.26, Rz 5.57, i?4 6.79, Rw 13.46, etc... 


12 Continuity result 

Theorem ILX Fix p > N, (3 = 1, let "fQ{l3,p, H) = 70 > 0 be the constant given 
by Lemma 
and /i = 7 


2.3 


let 0 < 7 < min{7o, ^}, be the global minimizer of (1.12) 
Note that w : pcfl —>■ K given by w{x) := Uj{p,~^x) is a weak 
solution in 77 of — pf'/S.w = w — w^ in pFl. Also note that p > '/S if 

(n) < 1 


4.1 


and Lemma 2.3 


we have that IIm.^ 


< C for some C > 6 independent of 7. 


lL“(n) = \Ml- 


By Proposition 
and ||M7l|c5.°(n) 

Let u* € 77q(H) be a global minimizer of (1.12) in 77g(f2) with 7 = 0. It is 


well known that u* is a unique global minimizer, smooth, strictly stable, and it 
does not change sign (see, for example, 0). 

Now, since is bounded in independently of 7 it is easy to see that 
u-y —>■ u* in as 7 —>■ 0, by the uniqueness of global minimizers of the limit 


problem (1.11) with 7 = 0. 


Let G € C^(K X C''^(H)) be given by G{y,u) = — Au — u + u^. Notice 

that duG{"f,u) G £(C'"^(H),K) and duG{0,u*) has trivial kernel by the strict 
stability of u*. Therefore, by the implicit function theorem (see for example 
[5U1 Theorem 1.1.1]), there is a neighborhood / x F C M x C"*(r2) of (0, u*) and 
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a continuous function X : I ^ V with A(0) = u* such that G{j, A(7)) = 0 for all 
7 € / and every solution of G{j, m) = 0 in / x 1/ is of the form (7, A(7)) for some 


'y G I. Since 


in G'^ as 7 —>■ O’*" and Uj is an arbitrary global minimizer, 


we obtain that is the unique global minimizer for all 7 € /. Finally, if the 


first Dirichlet eigenvalue Ai(n) < 1, then u* ^ 0 (cf. proof of Theorem 5.2) and 


by the strong maximum principle u* > 0 in fl, by the Hopf Lemma d^u* < 0 in 
(9n, where ly denotes the exterior normal vector, and therefore > 0 in 17 for 
all 7 G /, by making I a smaller neighborhood of 0 if necessary. 


□ 

Remark 12.1. Note that the proof of Theorem o also shows the existence of 
solutions in C'^’“(17) for equation (1.11) with 7 G [—70,0]. 


13 Numerical approximations 

In this Section we present some numerical approximations of a bifurcation 
branch with respect to jS from ft = ^ for the problem ( |1.1[ ) (cf. Theo¬ 

rem 11.11. These approximations were computed with the software AUTO-07P 



Figure 4: Here = (0, 2n). In this picture we present approximations of solutions of 
along the branch. Note that for ^ < 0 the positivity is lost and even negative solutions appear. 
Here the branch does not return to the region /3 > 0. 
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Figure 5: Here n = (0, IOtt). In this picture we show numerical evidence that positive 
solutions are not bounded by one and present oscillations. Also that the bifurcation branch 
may return to positive values of /3 (alth ough the positivity is lost). On the right we have an 
approximation of two solutions of along the branch corresponding to ^ = 0.35. 
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